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【摘要 】 风 荷载 往往 是 超 高 层 建筑 结构 安全 性 和 舒适 性 的 控制 性 荷载 。 本 文通 过 大 涡 模 拟 技术 并 结合 作者 建议 
的 一 种 新 的 满 流入 口 生 成 方法 (NSRFG 方法 ) ,对 四 种 不 同 锥 度 的 模 形 建筑 模型 进行 了 风 效 应 的 数值 模拟 对 比 研 
究 , 以 比较 不 同 锥 度 对 结构 风 荷 载 和 风 振 响应 的 影响 ,及 检验 本 文 数值 模拟 方法 的 适用 性 。 数 值 模拟 结果 表明 , 虽 
然 不 同 锥 度 的 建筑 模型 平均 风 压 和 脉动 风 压 分 布 规律 相似 ,但 在 结构 响应 方面 ,结构 横 风 向 峰值 基底 弯 矩 响应 随 
VO 模型 锥 度 增 大 显著 减少 ,表明 采用 适当 锥 度 的 体形 可 以 有 效 地 减 小 结构 横 风向 气动 荷载 。 本 文 数 值 模拟 结果 和 风 
过 洞 试验 结果 整体 规律 相符 ,表明 结合 改进 的 入 流 涡 流 生成 技术 的 大 涡 模 拟 具有 一 定 的 精度 ,可 以 给 超 高 层 建筑 气 


< 动 外 形 优化 研究 提供 有 价值 的 参考 。 
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1 二 前 
忆 3 随 着 社会 经 济 的 发 展 , 国 内 近年 来 已 出 现 了 越 
来 工 多 的 超 高 层 建筑 。 由 于 高 层 建筑 结构 具有 柔 
诬 类 .质量 轻 、 自 振 频 率 以 及 阻尼 较 低 的 特点 ,结构 
的 由 效应 逐步 成 为 控制 高 层 建筑 安全 性 、 舒 适 性 和 
经 济 性 的 最 重要 的 因素 之 一 。 超 高 层 建 筑 结构 属 
风 敏 感性 结构 ,对 风 荷 载 所 造成 的 扰动 十 分 敏感 ， 


随 着 计算 风 工 程 在 国内 外 的 发 展 ,数值 模拟 方 
法 已 经 成 为 结构 风 工 程 的 主要 研究 方法 之 一 ,大 涡 
模拟 (Large Eddy Simulation, LES ) 在 建筑 风 效 应 的 
模拟 中 得 到 越 来 越 广泛 的 应 用 。 随 着 数值 模拟 技 
术 和 计算 机 性 能 的 提升 ,近年 来 研究 人 员 尝 试 将 大 
涡 模 拟 运用 到 工程 实践 中 。 例 如 , 李 秋 胜 等 人 
采用 LES 分 析 了 高 层 建筑 的 空气 动力 学 效应 ,Dag- 
new, À. 等 人 0] 和 Daniels, S. J. 等 人 [4 运用 LES 对 
高 层 建筑 的 风 荷 载 进行 了 评估 。 


在 强风 作用 下 容易 受到 损坏 。 强 风 作 用 下 超 高 层 
建筑 结构 的 风 荷载 是 结构 安全 性 和 和 舒适 性 的 控制 
荷载 ,因此 对 风 荷 载 的 研究 和 分 析 已 成 为 超 高 层 建 
筑 结构 设计 和 体型 优化 的 重要 环节 。 

建筑 外 形 是 决定 作用 于 结构 上 气动 荷载 的 关 
键 因 素 , 通 过 合理 改变 建筑 物 的 外 形 , 对 建筑 外 形 
进行 优化 设计 ,可 以 显著 改变 作用 于 建筑 结构 的 风 
荷载 和 风 致 响应 。 通 过 传统 的 增强 结构 的 强度 与 
刚度 已 经 很 难 满足 现今 结构 设计 者 对 于 抗 风 设计 
的 要 求 ,于 是 合理 地 优化 作用 于 超 高 层 建筑 结构 上 
的 风 荷 载 已 成 为 结构 风 工 程 的 一 个 重要 研究 方向 ， 
如 台北 的 101 大 厦 及 阿联酋 迪拜 大 厦 的 外 形 都 进行 
了 抗 风 优化 研究 。 


【作者 简介 】 


目前 ,采用 LES 数值 模拟 方法 对 超 高 层 建筑 的 
气动 外 形 进 行 参数 优化 研究 还 较 少 ,分 析 可 能 原 
,由 于 超 高 层 建 筑 风 荷载 效应 对 气动 外 形变 化 较 
敏感 ,这 需要 数值 模拟 具有 很 高 的 精度 和 分 辨 率 ， 
而 通常 数值 风 洞 技术 还 很 难 做 到 和 物理 风 洞 类 比 。 

本 文 旨 在 运用 计算 流体 动力 学 (Computational 
fluid dynamics , CFD) 数值 模拟 技术 ,结合 作者 建议 
的 一 类 新 的 入 流 漠 流 生成 技术 ,对 不 同 棉 形 外 形 的 
超 高 层 建 筑 模 型 进行 大 涡 模 拟人 研究 ,以 评估 不 同 锥 
度 横 形 外 形 建筑 模型 的 气动 荷载 和 风 致 响应 影响 规 
律 ,并 探讨 采用 新 的 和 人流 清流 生成 方法 的 LES 数值 模 
拟 技术 研究 超 高 层 建筑 气动 外 形 优化 的 可 行 性 。 


余 远 林 (1991 - ) , 男 ,硕士 ,主要 从 事 建筑 结构 抗 风 和 数值 模拟 研究 。 
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2 数值 模拟 概况 


2.1 算 例 模型 与 结构 动力 特性 
本 文选 取 4 个 不 同 外 形 的 建筑 模型 进行 计算 ， 

模型 高 度 均 为 910mm , 模型 底 边 宽度 为 100mm , 缩 
尺 比 取 为 1:600 ,相当 于 实际 建筑 高 度 546m ,建筑 
高 宽 比 9.1。4 个 模型 中 包含 1 个 正方 形 等 截面 模 
型 (模型 1) , 3 个 锥 度 (模型 底部 截面 边 长 与 顶部 截 
面 边 长 差 值 与 模型 高 度 的 比值 ) 分 别 为 2.2%， 
4. 496 ,6.6 纺 的 枫 形 建筑 模型 (模型 2 ,模型 3 ,模型 
4) 各 模型 形状 和 尺寸 如 图 1 所 示 。 为 方便 叙述 tE 
度 用 表示 并 定义 : 

pam car (1) 

eo H 

stb iD, 和 D, 分 别 为 模型 底部 和 顶部 截面 边 
KSH 为 模型 高 度 ; 


ER 


Ln ) 


模型 1 模型 2 模型 3 模型 4 


图 1 模式 形状 和 尺寸 (mm) 


建筑 模型 动力 计算 所 采用 的 前 两 阶 自 振 周 期 
均 为 10. 2s( 自 振 频 率 0. 098 Hz) ,假定 振 型 沿 高 度 
方向 线性 分 布 ,阻尼 比 取 为 3.5% ,基本 风 压 参照 深 
圳 地 区 ,选取 重 现 期 100 年 .50 年 和 10 年 的 基本 风 
压 分 别 为 0. 90kPa .0. 75kPa 和 0. 45kPa。 

2.2 灌流 入 口 的 生成 

入 口 满 流 的 模拟 一 直 以 来 是 大 涡 模 拟 的 关键 

问题 之 一 ,其 意义 等 同 于 物理 风 洞 试验 中 需 首 先 准 


确 模 拟 边 界 层 风 场 特性 ,因此 有 许多 专门 研究 针对 
入 口 涪 流 的 准确 模拟 、 以 及 入 口 沸 流 影响 等 进行 深 
入 探 讨 。 本 文采 用 作者 建议 的 一 种 潮流 入 口 生 成 
方法 一 一 NSRFG 方法 中 生成 大 涡 模 拟 的 油 流 入 口 ， 
其 表达 式 如 下 式 所 示 : 

ux.) = Y p,sin(k, "T+2nft + 8,) (2) 

其 中 ,wu = {u,v,w},X = ru = x,y,z}; p, 
J| 2SCP. ) Afs S = fs ,5 5, 1 为 三 个 方向 上 的 风速 谱 


CRER) ,N 为 功率 谱 离散 数目 , 凡 = P Arr 
为 带宽 , o,- UC- m,r) ,另外 参数 取 值 由 下 式 得 


到 : 


U 
L, = m3 (3) 
yGf, 
-0.3 
其 中 ,yy = es B < n B » CD/Lu, C 为 
2.1,8 2 6.0 


衰减 系数 ; Lu X u Zr p iti D £A) NJ; D 20.5 ~ 
1. Oh, h 为 建筑 高 度 ; Uav 为 平均 风速 。 

通过 编程 将 NSRFG 方法 写成 UDF 用 户 自 定义 
程序 并 与 商用 软件 包 ANSYS FLUENT 链接 进行 
计算 。 

本 文 研究 对 象 原型 为 500m 以 上 的 超 高 层 建 
筑 , 由 于 此 高 度 已 经 超过 了 我 国 现行 荷载 规范 的 梯 
度 风 高 度 , 故 没有 采用 《建筑 结构 荷载 规范 
(GB50009-2001)》 的 风 场 参数 取 值 ,而 是 依据 英国 
工程 科学 数据 库 ( ESDU ) 的 方法 定义 了 数值 模拟 入 
口 的 平均 风 前 面 和 淇 流 度 剖面 ,如 图 2 所 示 。 
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2.3 网 格 划分 与 参数 设 定 
数值 模拟 在 华南 理工 大 学 风 洞 实验 室 的 数值 
风 洞 平台 商用 CFD 软件 包 ANSYSFLUENT 上 进行 。 
参照 AJ "和 相关 文献 的 规定 ,确定 流域 尺寸 如 图 
3a 所 示 ,堵塞 率 为 0.8% ,满足 小 于 3% 的 要 求 。 边 
条 件 方面 ,模型 表面 和 地 面 设 置 为 无 滑 移 壁 面 ， 
流域 左右 壁面 和 顶 面 设置 为 对 称 面 ,流域 人 口 为 速 
度 入 口 ,出 口 采用 压力 出 口 ,如 图 3a 所 示 。 由 于 建 
筑 模型 较为 规整 ,本 次 数值 模拟 采用 结构 网 格 进行 
计算 。 考 虑 到 计算 精度 和 计算 资源 的 平衡 ,最 小 网 
格 尺 寸 取 为 0.06D(D 为 模型 底面 宽度 ) ,网 格拉 伸 
比例 为 1.05。 近 墙壁 y+ 值 约 为 30 ~75, 适 用 于 壁 
Hi REL 不 同 模型 网 格 划分 方式 一 致 ,网 格 数量 约 
-— 万 ,网 格 划分 情况 如 图 3b 所 示 ( 仪 给 出 模型 
青 况 ) 。 
A ISO ASUES 隐 式 分 割 算法 
( Piéssure Implicit with Splitting of Operators, PISO) , 
离散 采用 中 心 差分 格式 ,时 间 离 散 采 用 二 阶 隐 
式 。 亚 格子 模型 采用 壁面 自 适 应 局 部 涡 粘 模 
型 Wall-Adapting Local 
ALE) 。 计 算 在 华南 理工 大 学 风 洞 实验 室 的 64 核 
号: 工作 站 进行 ,采用 16 核 并 行 计 算 。 时 间 步 长 
XA 0. 002s ,数值 模拟 计算 总 时 间 108 ,并 取 后 8s 的 
作为 分 析 数 据 。 参 考 高 度 采用 模型 高 度 
m ,在 参考 高 度 处 的 风速 为 11m/s。 人 参考 高 度 
处 建筑 模型 雷诺 数 约 为 7. 8 x 10* 。 


Eddy-viscosity model, 


(a) 计算 流域 尺寸 
图 3 计算 流域 尺寸 与 网 格 
3 结果 分 析 


Dp a 
脉动 风 压 系数 分 布 , 如 图 5 所 示 。 由 图 5 可 得 ,拥有 
不 同 锥 度 的 模型 平均 风 压 系数 和 脉动 风 压 系数 大 


(b) 网 格 划 分 


小 差别 不 大 ,分 布 相似 。 随 着 锥 度 的 增 大 , 风 压 分 
布 图 形 趋 于 拉 长 。 值 得 注意 的 是 ,在 模型 背风 面 肪 
动 风 压 系数 分 布 方面 , 随 着 锥 度 的 增 大 ,图 形 辐射 
中 心 由 两 个 化 为 一 个 , 且 有 上 移 的 趋势 。 具 体 风 压 
系数 分 布 情况 见 图 5。 

图 6 给 出 模型 1 在 10s 时 刻 满 流 场 示 意图 。 由 
图 6 可 见 ,来 流风 在 建筑 迎风 两 侧 发 生 分 离 ,并 产生 
分 流 涡 旋 , 在 来 流风 的 作用 下 交替 向 下 游 汇 出 , 产 
生 规 则 的 洲 涡 脱落 现象 。 模 型 的 尾 流 流 场 出 现 有 
别 于 来 流 滑 流 的 复杂 特征 注 流 ,呈现 出 典型 的 钝 体 
分 离 流动 特性 。 
3.2. 横 风 向 基底 气动 弯 矩 和 结构 响应 

强风 作用 下 风 和 荷载 将 成 为 影响 结构 安全 的 控 
制 荷 载 , 超 高 层 建筑 的 横 风 向 风 致 振动 响应 通常 会 
超过 顺风 向 响应 而 成 为 结构 设计 中 要 考虑 的 
EA, 


(pb) 背 风 面 平均 风 压 系数 
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(d) 迎 风 面 脉动 风 压 系数 
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1 f Ul 
模型 2 模型 3 


S (e) 背 风 面 脉动 风 压 系数 (人 侧面 脉动 风 压 系数 
o 5 风 压 系数 分 布 


q 


01234567 8 9101112131415 


(a) 中 垂 面 速度 云图 
6 流 场 : 10s 时 速度 云图 和 流 线 ( 模 型 1) 


1234567 8 91011121314 


(b) 1/2H 水 平 切面 速度 云图 


当 结 构 的 振 型 沿 高 度 线性 分 布 时 ,基底 弯 和 矩 和 
结构 的 1 阶 广义 力 成 正比 ,因此 基底 弯 矩 是 衡量 结 
构 整 体 气动 性 能 的 重要 指标 。 为 研究 不 同 锥 度 直 对 
模型 基底 弯 矩 的 影响 ,在 数值 模拟 过 程 中 记录 每 一 
个 时 间 步 的 横 风 向 基底 弯 矩 系数 Cy ,其 定义 如 式 
(4) Brzn : 

M 


Ci eu (4) 
30V, DH 


其 中 ,Ma 为 建筑 模型 的 横 风 向 基底 弯 和 矩 , 太 为 


参考 高 度 风速 ,D 为 建筑 模型 的 宽度 (此 处 取 底 部 
宽度 ) ,五 为 建筑 模型 高 度 。 


图 6 给 出 了 不 同 锥 度 的 模型 横 风 向 基底 弯 抑 
系数 功率 谱 密度 的 对 比 。 由 图 6 可 见 ,气动 弯 矩 随 
模型 截面 的 收缩 而 衰减 ,同时 基底 弯 矩 的 主峰 频率 
PESEE E 的 增 大 而 增 大 。 此 结果 与 文献 17] 的 风 洞 
试验 结果 基本 一 致 ,进一步 验证 了 本 文 数值 模拟 方 
法 的 适用 性 和 准确 性 。 

Across Wind Moment Coefficient C,,. 
— ER 
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图 6 不 同 锥 度 下 模型 横 风 向 基底 弯 矩 功率 谱 密度 


随 着 & 的 增 大 ,基底 弯 算 的 主峰 频率 也 增 大 ,使 
得 建筑 的 基底 气动 弯 矩 的 主 能 量 分 布 向 高 频 端 移 
动 , 这 可 能 会 对 结构 带 来 不 利 的 影响 。 图 6 中 竖 线 
对 应 的 频率 值 分 别 为 100 年 10 年 重 现 期 下 原型 结 
构 1 阶 自 振 频 率 换 算 到 模型 的 值 ,其 换算 关系 为 : 

f fA 

LE u^ (5) 

其 中 ,f, MS, 为 模型 与 原型 频率 ; AAA 分 
别 为 原型 和 模型 的 频率 比 缩 尺 比 和 风速 比 。 

根据 结构 动力 学 随机 振动 理论 ,由 式 (4) 记录 
的 基底 弯 和 矩 时 程 和 结构 动力 特性 可 计算 得 到 考虑 
结构 动力 放大 效应 的 横 风 向 基底 弯 矩 功率 谱 密 度 
和 标准 差 ,计算 公式 如 式 (6) ~ (8) 所 示 : 


Sy,() «1 HG) A (6) 
ow = | Sm df ri + nS ID 
(1) 
Ino r = 
Ur -G- 4T E48 QI 
(8) 


式 中 , | HO) V 为 结构 机 械 导 纳 函 数 ; fs LC 
分 别 为 结构 1 阶 自 振 频 率 和 阻尼 比 。 式 (7) 中 的 
ow FIL afe Su Cf )/ (AC) 分 别 为 弯 矩 响应 的 背景 分 量 
和 共振 分 量 。 在 得 到 基底 弯 和 矩 标准 差 之 后 ,模型 的 
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M max = M tg . O yp (9) 


式 中 , M, 为 峰值 基底 弯 矩 响应 ; M 为 基底 弯 
和 抢 响 应 平均 值 ( 横 风 向 取 为 0) ; g 为 峰值 因子 ,由 下 
式 可 得 : 

g, = V/2ln(v,T) +0.577/ J/2ln(v,T) (10) 

其 中 ,T 为 统计 时 间 ,通常 取 6008; v, 是 单位 时 
间 内 平均 位 移 的 水 平 交 叉 率 数 ,由 下 式 求 得 : 


v, = [JP S, dn/ [S4 dn]? (11) 


关于 高 层 建筑 结构 峰值 响应 计算 的 具体 方法 
详 见 文献 [8] 。 本 文 计算 响应 用 的 结构 自 振 周期 是 
根据 已 有 工程 和 规范 建议 估计 的 ,但 实际 的 建筑 结 
格 周 期 会 在 这 个 基础 上 发 生 一 定 范围 内 变化 。 为 
全 页 考虑 建筑 外 形 对 结构 气动 性 能 的 影响 必须 分 
枚 该 结构 风 致 响应 对 周期 变化 的 敏感 性 。 图 7 给 出 
了 3R 同 & 值 下 模型 峰值 基底 弯 矩 随 周期 比 的 变化 曲 
线 疼 。 此 处 周期 比 的 定义 是 实际 建筑 结构 周期 与 
布 文 假定 的 计算 周期 之 比 ,周期 比 的 变化 范围 为 
0:67- 1.4, 相 当 于 结构 的 计算 周期 范围 为 6. 12 ~ 
1 人 28s, 这 一 范围 基本 覆盖 了 《建筑 结构 荷载 规范 
(GB50009-2001) 》 计 算 结构 周期 的 经 验 公式 范围 ， 
基 太 满足 此 类 高 度 超 高 层 建 筑 自 振 周期 的 要 求 。 
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从 图 7 可 得 , 当 周 期 比 大 于 等 于 0. 85 时 ,峰值 
AUI 2S 4B BG HERE E 的 值 增 大 而 减 小 ,由 此 看 来 建 
筑 模型 的 截面 收缩 可 十 分 有 效 地 降低 结构 峰值 基 
底 弯 和 矩 。 当 周期 比 小 于 0.7 时 ,上 22.296 88 E = 
4. 496 两 种 外 形 模型 的 峰值 基底 弯 徐 反而 出 现 超过 


等 截面 模型 的 情况 ,可 见 变 化 规律 与 图 6 所 示 基 底 
弯 和 矩 功率 谱 密度 的 变化 规律 一 致 , 且 在 所 关心 的 周 
期 比 变化 范围 内 ,& = 6.6% 模型 的 峰值 基底 弯 矩 响 
应 最 小 。 此 计算 结果 与 文献 [8] 风 洞 试验 结果 变化 
规律 整体 相符 ,证 明 结 合 改进 的 入流 汕 流 生成 技术 
的 大 涡 模 拟 具 有 相当 的 精度 ,可 以 给 超 高 层 建筑 气 
动 外 形 的 优化 研究 提供 有 价值 的 参考 ,对 不 同 外 形 
的 建筑 模型 风 效 应 进行 定性 的 比较 和 判断 。 

但 不 可 否认 的 是 ,数值 模拟 结果 与 试验 对 比 仍 
存在 一 定 差距 ,可 能 与 大 涡 模 拟 的 滤波 效应 数值 
模拟 参数 设置 等 因素 有 关 , 将 在 后 续 的 工作 中 进行 
详细 比较 分 析 , 以 进一步 提高 数值 模拟 结果 的 精度 。 


4 结论 


本 文采 用 CFD 大 涡 模拟 技术 并 采用 作者 建议 
的 一 种 湛 流 生成 方法 一 一 NSRFG 方法 对 拥有 不 同 
锥 度 的 横 形 外 形 超 高 层 建 筑 模型 进行 了 风 效 应 的 
研究 ,得 出 以 下 结论 : 

(1) 不 同 模型 的 平均 风 压 系数 和 脉动 风 压 系数 
大 小 差别 不 大 ,分布 规律 相似 ; 随 着 锥 度 的 增 大 , 风 
压 分 布 图 形 趋 于 拉 长 。 

(2) 结 构 横 风向 峰值 基底 弯 矩 响应 随 模 型 锥 度 
增 大 显著 减少 ,这 说 明 超 高 层 建 筑 采 用 由 下 到 上 逐 
渐 收 缩 的 外 形 可 以 有 效 地 减 小 结构 横 风 向 气动 荷 
载 。 一 定 周期 比 范围 内 , 锥 度 越 大 , 横 风 向 气动 弯 
和 矩 的 削减 效果 越 好 。 锥 度 分 别 为 2.2% ,4.4%， 
6. 696 的 建筑 模型 在 周期 比 为 1 的 情况 下 ,相对 于 等 
截面 模型 横 风 向 基底 弯 矩 响应 减少 了 58.896, 
84. 5% 和 86. 296 。 同 时 锥 度 越 大 ,建筑 的 洲 涡 脱落 
频率 也 会 相应 增 大 ,这 有 可 能 造成 在 结构 自 振 频 率 
高 或 者 低 风 速 情 况 下 , 锥 度 高 的 建筑 风 致 响应 比 等 
截面 还 大 的 情况 。 此 类 建筑 的 外 形 优化 还 需 和 和 角 
区 优化 .避难 层 开 透 风 槽 等 局 部 气动 优化 措施 综合 
考虑 。 

(3) 采 用 作者 建议 的 新 的 LES 入 流 满 流 生成 方 
法 NSRFG 方法 ,数值 模拟 结果 与 风 洞 试验 客观 上 虽 
存在 一 定 差距 (可 能 与 大 涡 模拟 滤波 效应 和 数值 模 
拟 参 数 设置 有 关 ) ,但 仍 具 有 相当 的 精度 和 分 辩 率 ， 
可 以 较为 准确 地 预测 模型 流 涡 脱落 频率 ,给 超 高 层 
建筑 气动 外 形 的 优化 研究 提供 有 价值 的 参考 ,可 以 
达到 在 初步 设计 阶段 对 不 同 外 形 的 建筑 模型 风 效 
应 进行 体型 优选 的 目的 。 
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LES Research on Wind Effect on Tapered Super-tall Buildings 


Yu Yuanlin, Yang Yi, Xie Zhuangning 


(State Key Laboratory of Subtropical Building Science, South China University of Technology , 
Guangzhou 510641, China) 

C Abstract ; Wind load was often the key load for the safety and comfort of super-tall building structures. Large 
eddy simulation ( LES) technique combined with a new inflow turbulence generator( NSRFG) was employed to study 
the wind effects on tapered building models with different tapers, which aimed to compare the effect of different ta- 
per on wind loads and wind-induced responses of the structure and validate the applicability of the numerical meth- 
od of this paper. Numerical simulation showed thatthe distributions of the average and fluctuating pressure on build- 
ing surfaces were similar, whilefor the wind-induced response, the peak base bending moment ( BBM ) of different 
models in across wind direction decreased remarkably with the increasing of the taper. Therefore, a building with 
bottom-up gradually shrinking and an appropriated tapershape can effectively decrease the aerodynamic wind load in 
across-wind direction. In addition , results from the numerical simulation were coincided with that from the wind tun- 
nel test in the whole , which indicated that LES technique combined with an appropriated inflow turbulence generator 
would havecertain accuracy and canprovide valuable reference to the aerodynamic shape optimization of super-tall 
buildings. 

Key Words :High-rise Building; Aerodynamic Optimization; Large Eddy Simulations; Wind Loads; Structural 
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